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ABSTRACT 

We present a method for investigating variations in the upper end of the stellar Initial Mass Function 
(IMF) by probing the production rate of ionizing photons in unresolved, compact star clusters with 
ages < 10^ yr and with different masses. We test this method by performing a pilot study on the young 
cluster population in the nearby galaxy NGC5194 (M51a), for which multi-wavelength observations 
from the Hubble Space Telescope are available. Our results indicate that the proposed method can 
probe the upper end of the IMF in galaxies located out to at least ~10 Mpc, i.e., a factor w200 
further away than possible by counting individual stars in young compact clusters. Our results for 
NGC5194 show no obvious dependence of the upper mass end of the IMF on the mass of the star 
cluster down to wlO'^ M©, although more extensive analyses involving lower mass clusters and other 
galaxies are needed to confirm this conclusion. 

Subject headings: galaxies: star clusters: general - galaxies: star formation - galaxies: individual 
(NGC5194) - stars: massive - stars: mass function 



1. INTRODUCTION 

The distribution of stellar mass es at birth (th e stel- 
lar Initial Mass Function, or IMF. ISalpeteii HOSSl ) is the 
fundamental property, together with the efhciency of star 
formation, that quantifies the conversion of gas into stars 
in galaxies. It enters, directly or indirectly, in the deter- 
mination of many physical quantities of stellar popula- 
tions and galaxies, such as the conversion of luminosity 
to mass and star formation rate (SFR). The IMF is often 
parametrized by the broken pow er law form of [Kroupal 
I2OO2) or the log-normal form of lChabrieJ (j2003l i^ 

Although many studies indicate that the stellar IMF 
is consistent with a universal form (see, for a review, 
iBastian et al.l [2010), a few recent investigations cast 
doubts on the universality in all environments. 

By comparing galaxy-integrated properties of large 
samples of nearb y gala x ies at multiple wave l engths , 
Hovcrsten & Glaze brookl (I2008D. iMeurer et all (|2009f ). 
Lee et all ([2009), and'BoseUi et al.' (2009) find that, on 



average, the ionizing/non-ionizing flux ratio of galax- 
ies decreases for decre asing gala. x y luminosity or mass 
or SFR or SFR/area. iLee et all (|200al exclude, based 
on statistical arguments, that galaxy-integrated stochas- 
tic sampling or variations in the star formation histories 
alone can fully account for the observed trend. A sys- 
tematic variation of the high end of the IMF with some 
galactic or local parameter (luminosity, SFR, SFR/area, 
gas density) has been suggested as one of the possible 
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ways to explain th e obs ervations ([Krumholz fc McKe 
20081: IMeurer et al.ll2009l : lPfiamm-Altenburg et al.ii2009 
Lee et al.l I2009D. althoii g h oth er scenarios ar e pos 
sible (e.g., [D ong ct a ll l2008t 'Kotull a et all 120081: 



Gogarten et al.i 2009: Hunter et al. 2010: Goddard et all 
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The model suggested by iPflamm-Altenburg et al.l 
(|2009f l invokes the existence of two correlations: one be- 
tween a galaxy's SFR and its clusters' maximum mass 
(e.g.,|Larsen 2002), and another (deterministic) between 
a cluster mass an d the most massive possi ble star the 
cluster can form (jWeidner fc Kroupal 120061) . The na- 
ture of the second correlation implies that clusters are 
not subje ct to th e effect of stochastic sampling of the 
IMF (El megreenl [2006), i.e. 100 clusters each with 
mass=10^ M0 would not be equivalent to one 10^ M© 
cluster, as the former will systematically lack the high- 
mass stars contained in the latter. 

The standard approach to deriving the functional form 
of an IMF, i.e. counting the individual stars con- 
tained in single-age, young stellar populations (typically 
<3 Myr old star clusters), is difficult to apply over sta- 
tistically large samples of clusters at a given mass and in 
a variety of galactic environments, and provides limited 
leverage for discriminating stochastic sampling from sys- 
tematic trends in the IMF. Significant blending will occur 
for massive stars due to crowding in the centers of star 
clusters (|Maiz-Apellanizl2008HAscenso et al.l2009i) . even 
with the angular resolution of the Hubble Space Tele- 
scope, for distances larger than those of the Magellanic 
Clouds (~55 kpc). These difficulties may acco unt for the 
contra sting results obtained by Maschberge r fc Clarkd 
(|2008f ). who concluded that the most massive stars 
in Milky- Way clust ers are stochastically drawn from a 
universal IMF, and iWeidner. Kroupa fc Bonnelll ([2010) . 
who claimed evidence for a systematic trend between the 
masses of clusters and those of their most massive stars. 

Here we present an alternate method that exploits the 
integral properties of star clusters to place limits on vari- 
ations of the high-end IMF. Each star cluster is treated 
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as single unit, and IMF variations are investigated by 
measuring the dependence of the ionizing photon rate, 
normahzed by the cluster mass, on cluster mass (sec - 
tion 2). This method extends the iCorbelli et alJ (|2009l ) 
approach (the cluster birthline) by normalizing the ion- 
izing photon rate to the age-independent cluster mass, 
instead of the age-dependent bolometric luminosity. Our 
method, however, requires that accurate ages for each 
cluster be determined. We test our approach using obser- 
vations of the young star cluster population in the nearby 
galaxy NGC5194 (section 3). A summary of findings and 
future directions are given in section 4. 

2. METHOD AND PREDICTIONS 

Stellar population synthesis models provide specific 
predictions for the ionizing photon rate, N(H°), produced 
by a single-age stellar population as a function of time, 
for a given IMF, mass, and metallicity. Compact star 
clusters represent a close approximation to the ideal- 
ized case of a single-age population. The rate N(H°) 
decreases with increasing age, as massive stars die out, 
and ionized gas generally becomes undetectable in clus- 
ters older than ~8-10 Myr (more than 100 times fainter 
than at 1-2 Myr age). Most HII re gions form ex pand- 
ing shells around the central cluster ()Martinlll998[ ). thus 
decreasing the surface brightness of the ionized gas emis- 
sion lines even further. For practical purposes, ionized 
gas emission will not be easily detected from star clusters 
older than w6-8 Myr. 

Measuring N(H°) from a star cluster at a given (young) 
age is equivalent to measuring the number of massive 
stars distributed according to an IMF. This, in models, 
includes fractional stars in the high end of the IMF, for 
low-mass clusters. Real clusters, however, do not contain 
fractions of stars, and, for masses<a fewxlO"' Mq, they 
will be subject to stochastic sampling of the IMF, even 
if the underly ing IMF is 'univ ersal', i.e., independent of 
cluster mass ([Elmegreenl 120061 ) . Stochastic sampling in- 
duces a dispersion on the number of stars at a given 
stellar mass that are present in a real cluster, and the 
dispersion increases for decreasing cluster mass. This, in 
turn, affects a number of observational parameters, in- 
cluding the iiiferred i nass and N(H°) from cluste rs (e.g., 
iCervifio et all [200l iCervifio fc Luridianal [200l . The 
scatter on the mass determination increases from 6% for 
a 10^ M© cluster to --50% for a 10^ M© cluster; the 
scatter on N(II°) increase s from ^10% to 80% for the 
same cluster mass range (jCerviho et al.l[2002[ ). In our 
approach, the impact of stochastic IMF sampling is min- 
imized by combining large numbers of same-mass, young 
clusters to a cumulative equivalent mass «10^ M©, since 
a 'universal' IMF would be fully sampled in a cluster 
of this mass, and stochastic uncertainties are reduced to 
10% or less in both M^; and N(H°). 

Figure[Tl left, shows the model expectations for N(H°) 
(as traced by the Ha luminosity, L(Ha)) per unit clus- 
ter mass, as a funct ion of cluster mass, for both cases 
of a 'universal' IMF (|Kroupall2002l ) with fractional stars 
and a cluster-mass-dependent IMF with the maximum 
possible stellar mass dep ending on the cluster mass 
(iWeidner fc Kroupal 120061) . Results similar to those of 
iWeidner fc Kroupal (|2OO60 would be obtained if the IMF 
became steeper at the high end for decreasing cluster 
mass. The L(IIa)/Mci ratio is averaged over the age 



ranges 1.5-5 Myr and 2-8 Myr for constant star forma- 
tion, in order to bracket the uncertainty in cluster ages 
found in most observational analyses. While our default 
model has metallicity 1.5 solar (to match NGC5194, see 
below), we show examples also at 1/3 and 3.5 solar metal- 
licity. The metal content of the ionized gas has, in fact, 
strong influence on the ratio L(Ha)/Mc;; this parame- 
ter will be critical when comparing star cluster popu- 
lations belonging to different galaxies, while its impact 
will be less important when comparing clusters at differ- 
ent masses within the same galaxy, under the reasonable 
assumption that no systematic trend exists between a 
cluster mass and its metal content. 

Stochastic sampling of the stellar IMF in low-mass 
clusters, combined with photometric uncertainties and 
reddening, can also affect age determinations from multi- 
wavelength data. Clusters younger than ~10 Myr 
that do not have ionized gas emission can be as- 
signed older ages; conversely, clus ters older than 10 Myr 
can be assigned ages ~ 6-9 Myr (jChandar et al.|[2"010bl : 
iFouesneau fc Lanqon 1 12010|) . Thus samples of young, 
low-mass clusters are 'polluted' by either addition of old 
clusters or loss of young members. The two may com- 
pensate each other to some degree, but this is still un- 
quantified. Furthermo re, masses may be underestim ated 
by up to a factor of 2 (jFouesneau fc Lanpon Il2010l) . We 
discuss both effects below, when applying our approach 
to actual data. 

3. APPLICATION: A TEST CASE 

We test our method on th e nearby s tar-forming 
galaxy NGC5194 (D~8.4 Mpc, PTdkl [19881) . for which 
extensive multi-wavelength, broad- and narrow-band, 
archival imaging from the Hubble Space Telescope ex- 
ists (GO-10452 and GO-10501). The galaxy repre- 
sents an ideal target for this test: it is nearly face-on 
(i '^20°), thus minimizing line-of-sight overlaps among 
clusters, and its c luster population has been exten- 
sively studied (e.g. iBastian et al.l l2005t iChandar et al.l 
l2010bH. The HST observations of NGC5194 include U 
(with WFPC2), BVI and Ua (with ACS) for the en- 
tire disk, and Pq;(A=1 8756 ^m, with NIC MOS) for the 
inner '---3 kpc radius (jScoville et al.l [200l[) . The galaxy 
has a mean metal abundance about twice solar, which 
remains super-solar a cross the entire disk up to R25 
(jMoustakas et al.ll2010t ). 

We use the sam ple of clusters identified by 
IChandar et al.l (j2010b[ ). and refer the reader to that pa- 
per for the sample selection and measurements. Descrip- 
tions of the selection and age-dating methods for star 
cluster populations in external galaxies can be found 
also in Whitmorc et al. (2010, applied to the Anten- 
nae Galaxies) and IChandar et al.l (j2010al applied to 
NGC5236). We summarize here some of the steps that 
are relevant to the present paper. For all clusters, 
ages and extinctions are determined from the x^ min- 
imization between the observed U,B,V,I and Ha (non- 
continuum-subtracted) photometry and predictions from 
the Br uzual fc Chariot ( 2003 , 2007 Update) stellar pop- 
ulation mo dels combined with a Milky- Way-type extinc- 
tion curve ([Fitzpatrickl 1 1 999l ) . The choice of the extinc- 
tion curve does not impact the results, because most 
curves are degenerate in the wavelength range from U 
to I (|Calzetti.,2001.) . Changing the population synthesis 
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models, e.g., to Starburst99 (jLeitherer et al.|[l999f ). also 
has minimal impact on the results in the case of young 
clusters (^10 Myr). The use of the non-continuum- 
subtracted Ha photometric point (which becomes a R- 
band point in the absence of ionized gas) adds leverage 
for separating young, extincted clusters from old, rela- 
tively unextincted ones, even in those cases where th e 
young clusters do not ionize gas (jChandar et al.|[2010bD . 
The typical uncertainty is a factor ~2 in the final age 
determination, and ~60% in the mass. The mass of 
each cluster is estimated from its extinction-corrected 
V band luminosity and the age-dependent mass-to-light 
ratio from the stellar population models, assuming a 
Kroupa/Chabrier stellar IMF. Changing the upper mass 
limit of the IMF from 120 to 30 M© changes the mass 
estimate by «2.5X, which is within the range of uncer- 
tainties we consider. We consider two mass bins, located 
at the two extremes of the sampled range in the clus- 
ter mass distribution: ~1.5xl0^ M© and ~4xl0* M©. 
They span about a factor 30 in mean mass and a factor 
300 in individual cluster masses. 

We inspect each cluster visually, and remove from the 
sample any cluster located in a sufficiently crowded re- 
gion to hamper Ha line measurements. Our final, visu- 
ally inspected, sample consists of 86 clusters in the mass 
range 0.5-2 10^ Mq and 57 in the range 2-15 10^ Mq, 
all younger than 10 Myr. Of the 86 clusters in the lowest 
mass bin, 39 (45%) are undetected in Ha down to our 
3 a limit of 7x10^^ erg s~^ (calculated in a 5 pixels ra- 
dius aperture = 10 pc). As mentioned in section 2, some 
young clusters can be assigned ages older than 10 Myr, 
and vice-versa. We address this issue by adopting two 
extreme scenarios: case (a): all clusters whose colors in- 
dicate young ages are included in the average of each 
mass bin; case (b): only clusters which have colors of 
young populations and also have Ha detections above 
3 a are included in the mass bin average. In order to ac- 
count for the possibility that cluster masses in the low- 
mass bin may be underestimated by as much as a factor 
2 (Foucsncau & Lancon 2010), we include clusters with 
masses as low as 500 Mq. Of the 57 clusters in the higher 
mass range, only 6 (11%) are below 3 cr in Ha. 

Photometric measurements of the hydrogen recombi- 
nation lines (continuum-subtracted Ha-|-[NII] and Pa) 
are performed at the positions of the clusters, with aper- 
ture sizes that scale according to the cluster's mass, 
i.e., according to the expected Stromgren radius of the 
HII region. Aperture sizes are chosen to be about 
0.35 Kstroemgrem becausc of crowdiug of the clusters, 
and age-dependent aperture corrections are implemented 
based on the few 'isolated' HII regions identified in our 
sample. The final photometry is visually inspected to 
check for contamination from neighboring sources; the 
local background, determined for each aperture from an- 
nuli 2 pixels wide around each region, is subtracted from 
each photometric measurement, in order to avoid con- 
tamination by diffuse emission not associated with the 
cluster. The uncertainty in the final Ha-|-[NII] fiux from 
uncertainties in the aperture correction range from about 
20% for the youngest and most luminous regions to about 
80% in clusters older than ~5-6 Myr. [Nil] contamina- 
tion in the narrow-band filter i s removed by using th e 
average galactic [Nil] /Ha ratio (iMoustakas et al.]l2010l ). 



The Pa image only covers a portion of the NGC5194 disk, 
and is used here to derive ionized gas extinction correc- 
tions region by region for the area of ov erlap between 
the A CS and NICMOS mosaics (see, also, iScoville et al.l 
|2001| ). for a total of 29 regions. The median values of the 
gas extinction corrections in each cluster mass bin are 
then adopted for the regions outside of the overlap area. 
This may lead to an overestimate of the dust extinction 
correction of Ha in each mass bin, as the mean extinc- 
tion shows a gra dient of decreasin g values towards the 
outer regions (Cal zetti et al.l [20051 ): this will not impact 
the relative comparison between mass bins, only their 
absolute values. 

In Figure [H right, the ratio L(Ha)/Mci=< 
L(Ha)>/<Mci> for each of the two cluster mass 
bins is shown for the two color-age cases (a) and (b), as 
cyan and magenta points, respectively. The error bars 
of the high-mass bin are dominated by photometric and 
aperture-correction uncertainties, while those of the 
low-mass bin contain a significant contribution from the 
stochastic uncertainty in cluster mass determinations. 
The absolute values of L(Ha)/Mci are similar for the 
two mass bins, within the 1 a uncertainty, for both sets 
of points. The absence of an obvious differential trend 
in L(Ha)/Mci for the two mass bins is consistent with 
a universal, stochastically-sampled, IMF. However, we 
cannot completely eliminate the scenario for which the 
maximum possible mass of a star increases with the 
cluster mass, as the data are only 1-2 a away from this 
model's expectations (triangles). 

The data follow the 2-8 Myr age average points (aster- 
isks) of the 1.5 solar metallicity stellar population mod- 
els, implying the following, non-exclusive, possibilities: 
(1) our criteria select clusters fairly uniformly distributed 
in the age range 2-8 Myr; (2) the models used in Fig- 
ure [l] right, have too low metallicity for NGC5194, and 
a higher metallicity model should be used (cf. Figure [U 
left panel); (3) significant, mass-independent, ionizing 
photon leakage is present in our regions, which we have 
not attempted to correct for in our analysis. 

4. CONCLUSIONS AND FUTURE WORK 

We have presented a method to test for variations of 
the high-end of the stellar IMF that exploits the inte- 
grated properties of star clusters, rather than resolved 
star counts. This approach can overcome statistical lim- 
itations in IMF determinations and average out effects of 
stochastic sampling by leveraging the large numbers of 
star clusters in galaxies. The ratio of the ionizing pho- 
ton rate (as measured by the extinction-corrected Ha 
emission) from young clusters to cluster mass provides a 
powerful tool to test for such variations. A basic require- 
ment for this method is to isolate the young (<10 Myr) 
star cluster population, and measure its ionized gas emis- 
sion. Accurate (~2X) ages can be obtained with high- 
angular-resolution (e.g., HST) multi-col or observations; 
a min imum set appears to be UBVRI (iChandar et al.l 
'2010b) , although the addition of the UV can add an ex- 
tinction leverage (Calzetti et al. 2005) to further reduce 
the age uncertainty. Ionized gas emission measurements 
also need to be corrected for dust extinction, requir- 
ing observations of two hydrogen recombination lines. 
Comparisons between galaxies further requires that the 
metallicity of each galaxy is known, as the ionizing pho- 
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ton rate is highly dependent on metaUicity. 

We have tested this method by performing a pilot 
study on the cluster population in the nearby galaxy 
NGC5194. Our analysis shows that, when a sufficiently 
large number of similar-mass clusters is summed to- 
gether to average out effects of stochastic sampling of 
the IMF, the mean Ha luminosity normalized by the 
mean cluster mass is relatively independent of the clus- 
ter mass, down to Md ^10'^ M©. This is consistent with 
an effectively universal IMF independent of the cluster 
mass. We do not find strong evidence that the maxi- 
mum possible mass of a star in a cluster decreases sys- 
tematically with decreasi ng cluster mass (as proposed by 
iWeidner fc Krou"pall2006D . However, this is not the most 
stringent test of the stochastic sampling model because 
clusters more massive than 10"^ Mq are expected to con- 



tain more than 1 0-type star. Clearly, statistics need to 
be expanded by including lower cluster masses, and clus- 
ter populations from other galaxies, spanning a range 
of internal and external environments. The proposed 
method enables tests of variations of the high-end IMF 
up to at least ~10 Mpc distance (NGC5194 is located 
at 8.4 Mpc distance), and possibly more, depending on 
the crowding of the cluster population. This distance 
range is at least a factor of 200 larger than what has 
been achieved so far with direct star counts. 



This research has made use of the NASA/IPAC Extra- 
galactic Database (NED) which is operated by the Jet 
Propulsion Laboratory, California Institute of Technol- 
ogy, under contract with the National Aeronautics and 
Space Administration. 
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Fig. 1. — LEFT: The expected ratio L(Ha)/M(,; as a function of the cluster mass M^,; both for a universal IMF (crosses and asterisks), 
and for a cluster— mass- dependent IMF (filled squares and triangles, Weidner & Kroupa 2006). The models are from Starburst99 (2007 
Update. ILeitherer et al.iil999i ). and assume averages over large numbers of clusters. The mean L(Ha)/M;,; ratios are averaged over two age 
ranges: 1.5—5 Myr (crosses and squares) and 2-8 Myr (asterisks and triangles), joined by vertical bars, under the assumption of constant star 
formation for each age range. We assume that clusters younger than 1.5-2 Myr are too dust-enshrouded to provide significant contribution 
at optical wavelengths. The black symbols are for 1.5 solar metallicity models; examples of higher (^3.5 solar, red) and lower (~l/3 solar, 
blue) metallicity are shown at a representative mass. RIGHT: The ratio <L(Ha)>/<M(,; > as a function of the mean cluster mass M^; 
in two mass bins for NGC5194, for clusters younger than ~8— 10 Myr (stars with 1 a error bars), including all clusters within each mass 
bin (cyan, bottom symbols), and including only clusters with Ho flux detections >3 cr (magenta, top symbols). The number of clusters 
in each mass bin is shown close to the symbols. Expectations from models with 1.5 solar metallicity and averaged over the 2-8 Myr age 
range are shown (asterisks and triangles), as they are the closest match to the NGC5194 high-mass data. 



